Spin-orbit interaction in a two-dimensional electron gas: a 577(2) formulation 
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Spin-orbit interaction is usefully classified as extrinsic or intrinsic depending on its origin: the 
potential due to random impurities (extrinsic), or the crystalline potential associated with the band 
or device structure (intrinsic). In this paper we will show how by using a SU(2) formulation the 
two sources of spin-orbit interaction may be described in an elegant and unified way. As a result 
we obtain a simple description of the interplay of the two types of spin-orbit interaction, and a 
physically transparent explanation of the vanishing of the d.c. spin Hall conductivity in a Rashba 
two-dimensional electron gas when spin relaxation is neglected, and its reinstatement when spin 
relaxation is allowed. Furthermore, we obtain an explicit formula for the transverse spin polarization 
created by an electric current, which generalizes the standard formula obtained by Edelstein and 
Aronov and Lyanda-Geller by including extrinsic spin-orbit interaction and spin relaxation. 
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I. INTRODUCTION: THREE SURPRISES 

Spin-orbit interaction in metals and semiconductors 
has attracted much interest in recent yearsi^. One rea- 
son for this is the possibility offered by such an interac- 
tion to achieve all-electrical control of the electron spin, 
as e.g. in the case of the spin Hall effect. The latter 
is a phenomenon in which a spin current appears in the 
direction transverse to a charge current driven through 
a sample. Our paper is motivated by several surprises 
and puzzles in the theory of the spin Hall effect in a two- 
dimensional electron gas (2DEG) with linear Rashba (or 
Dresselhaus) spin-orbit coupling (this will be referred to 
in what follows as the Rashba 2DEG). 

The first surprise was Sinova and co-workers' obser- 
vation that in the absence of disorder or interactions 
the spin Hall conductivity assumes the universal value 
<j sH = e/8ir 3 . This is in sharp contrast with the behav- 
ior of the regular electrical conductivity, which, of course, 
is infinite under the stated assumptions, due to the ab- 
sence of a scattering mechanism that may establish a 
steady state regime. The spin Hall conductivity - within 
certain limitations - depends neither on the strength, nor 
on the microscopic form of the spin-orbit coupling. 

The second surprise was that scattering from impuri- 
ties, while restoring a finite electrical conductivity, leads 
to another universal result, namely o~ sH — for arbi- 
trarily small impurity concentrations^—. This result can 
be traced back to a very specific form of the equation of 
motion for the spin polarization in the Rashba 2DEG, 



d, 



> t s y = -2majy , 



(1) 



where s y is the spin polarization in the y-direction, m 
the effective electron mass in the sample, a the Rashba 
spin-orbit coupling constant and j' the spin current po- 
larized along z and flowing along ?A~— . In a steady-state 
situation the time derivative of the spin polarization van- 
ishes and so does the spin current and the d.c. spin Hall 



conductivity. 

At this point a remark on the origin of spin-orbit cou- 
pling is in order. Spin-orbit interaction is usefully classi- 
fied as extrinsic or intrinsic depending on its origin, either 
in the potential due to random impurities (extrinsic), or 
in the crystalline potential associated with the band or 
device structure (intrinsic). In Refs. [ill and [l2l the spin 
Hall effect was calculated for a 2DEG ignoring intrinsic 
spin-orbit coupling but assuming an extrinsic one of the 
form 



f so = Axvy(x). P 



(2) 



where er = (a x , o~ y , o~ z ) is the vector of Pauli matrices, 
V(x) is the random impurity potential, and Ao the ef- 
fective Compton wavelength characterizing the extrinsic 
spin-orbit interaction strength. Two distinct contribu- 
tions to the spin Hall conductivity were identified, the 
skew scattering (ss) and side jump (sj) contributions: 



T sff 



a 8 /, with 
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(3) 



where n is the electron density, pp the Fermi momen- 
tum, I the mean free path, and vq the amplitude of the 
impurity scattering potential. 

The next logical step for the theory was to combine 
the two mechanisms, which of course are simultaneously 
present in experiments. This led to a third surprise^—. 
Namely the spin Hall conductivity was found to be inde- 
pendent of the Rashba coupling constant a, and yet its 
value was different from the one at a = 0. In other words, 
the spin Hall conductivity appeared to be a non-analytic 
function of the Rashba coupling constant a - its value 
exhibiting an unphysical jump from a — to arbitrarily 
small but finite a. 

The issue of the non-analytic behavior of the spin Hall 
conductivity for a — > was settled in Ref. [l6|. The "un- 
physical jump" was shown to result from the neglect of 
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spin relaxation - an effect that appears to second order in 
Aq and was therefore missed in the first-order treatments 
of Refs. IT3I - TT5I When spin relaxation was included^ the 
discontinuous jump of a sH at a = was replaced by 
a smooth evolution on the scale of a = — where 

PF\/t s t ' 

1/t s (1/r) is the spin (momentum) relaxation rate. 

However, an additional problem remains. The theoret- 
ical results of Refs. [l3l - fl6l agreed that the Rashba term 
suppresses the skew scattering contribution to the spin 
Hall conductivity for arbitrarily small a. However, no 
consensus has been achieved so far on the side jump con- 
tribution. While Refs. [H and [H claim a sIi = (l/2)crf , 
in Ref. 14 a vanishing spin Hall conductivity is predicted, 
a sH = 0. 

Here we reconsider the problem in order to settle the 
issue of the side jump. We will make use of the anal- 
ogy between spin-orbit coupling and a spin-dependent, 
i.e. SU(2), gauge fiel d 17 i 18 . We will start with a discus- 
sion of the equation of motion for the spin density which 
allows to determine the spin Hall conductivity in the 
presence of intrinsic and extrinsic mechanisms, with the 
result that a sIi = 0. Subsequently we will also present 
the theory based on a generalized Boltzmann equation, 
and will discuss in the appendices how this is microscop- 
ically derived from the Keldysh formalism. It will be 
shown how the simple dynamical argument for the van- 
ishing of the spin Hall current in the steady state, Eq. ((T|) , 
is generalized to the more complex case in which both 
intrinsic and extrinsic sources of spin-orbit coupling are 
present. 

Finally, it is worth mentioning that quantum and semi- 
classical kinetic equations have been widely used to dis- 
cuss spin-orbit coupling effects, see e.g. Il9i - [23l Our treat- 
ment, which we here specialize to the Rashba model, is 
actually valid for any system described by a Hamiltonian 
with linear-in-momentum (intrinsic) spin-orbit coupling, 
poss ibly slowly varying in space and time; as discussed in 
Ref. [23, a Zeeman field could also be included. As such, 
it is a treatment which in principle allows one to study 
the interplay between intrinsic and extrinsic spin-orbit 
coupling mechanisms in a wide range of systems. 



II. EQUATION OF MOTION FOR THE SPIN 
DENSITY 

In this section we will study the equation of motion for 
the spin density in a 2DEG including disorder, intrinsic 
and extrinsic spin orbit coupling with the Hamiltonian 

p 2 y2 

H= ^+a( Py o- x -p x o-y)+V(x)-^<rxVV(x)-p. (4) 

Here V(x) is the sum the potential created by the elec- 
tric field that drives a current through the system and 
a statistically fluctuating potential due to impurities, 
V"(x) = eE-x + 6V(x) . The fluctuations of the impurity 
potential are given by 5V(x)5V(x') = n^dfe — x'), n. 



being the impurity concentration and vq the scattering 
amplitude. In this paper we will assume weak spin-orbit 
coupling, in the sense that apF *C ep and {XqPf) 2 "C 1. 
Since both spin-orbit coupling terms are linear in the mo- 
mentum p they can be parameterized by a SU (2) vector 
potential, 

H=^ + V(,) + ±At Pi a a . (5) 

Here the summation over indices i and a is implied. Like 
the scalar potential V also the vector potential A* is the 
sum of a smooth and a fluctuating contribution, Af = 
A? + 5Af, with 

A? = -2mae laz - ^emXle a jiEj , (6) 

5A$ = -imAge^a^x). (7) 

The equation of motion for the spin density operator (op- 
erators will appear with a hat ") assumes the compact 
form 

dts a + dij? = e abc A\f t . (8) 
The spin density operator itself is defined as {K = 1 ) 

s a (x) = ±{a a ,S(x-r)}, (9) 

whereas jf are the components of the spin current oper- 
ator. In analogy to the conventional current operator we 
write them as a sum of a paramagnetic and a diamagnetic 
term, jf = [j P ]f + [j D ]f, with 

Op}" = ^ {Pi*a, S(x ~ f )} , $ D ]t = ^A^(x - f ) . 

The goal is to obtain the equation of motion for the 
(quantum and disorder) averaged spin density. While the 
l.h.s. of Eq. © is easily averaged, the r.h.s. involves the 
fluctuations of the vector potential and its correlations 
with the current density, 



d t s a + ddf = tabcA^jf + ^A h Sl . (11) 

We have not been able to determine the correlator on 
the r.h.s. of the equation exactly. An analysis carried 
out in the limits of weak spin-orbit coupling, weak disor- 
der, and slow spatial and temporal variations shows (see 
Appendix lAl for details) that 

eabJAfdjf = • (12) 

More precisely, the correlations of the vector potential 
with the paramagnetic current operator cancel with the 
correlations of the diamagnetic current operator: 



e abc 5A\[8 ]P ]1 = -e abc SAl[Sj D ]9 = -e abc A°\j s $ , (13) 
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where [j s j]f = Cicjcr^ Ej is the spin Hall current from the 
side jump mechanism. The averaged continuity equation 
is thus 

dts a + dijt = lahAii ■ (14) 

The argument for the vanishing spin Hall conductivity 
is thus the same as in the absence of extrinsic spin-orbit 
coupling: in a stationary and homogeneous situation the 
l.h.s. of the equation is zero. From the y-component of 
(fT4|) we obtain 2majy = 0, and thus the spin Hall con- 
ductivity is zero. 

We can gain a deeper insight into the physical meaning 
of this mathematical result by examining the structure 
of the Boltzmann equation in the presence of extrinsic 
and intrinsic SU(2) vector potentials. We will show in 
the following sections that the net spin Hall current due 
to an electric field in the x-direction is the sum of 

(i) A "drift" contribution, associated with the skew- 
scattering and side-jump due to impurity scatter- 
ing. 

(ii) A Hall current due to the "SU(2) magnetic field" 
(see Section 3.2) associated with the SU(2) vector 
potential of the Rashba spin-orbit coupling. 

(iii) A diffusion contribution, which is proportional to 
the y-component of the uniform spin polarization. 
The spin polarization originates from the so-called 
Edelstein effect, i.e., the fact that electrons drifting 
in the ir-direction in the presence of a Rashba field 
"see" an effective magnetic field pointing in the y- 
direction. 

The presence of a spin diffusion current arising from a 
spatially homogeneous spin distribution may appear puz- 
zling at first, but is a natural effect of the SU(2) vector 
potential. The point is that the spin distribution, while 
spatially uniform, is non-isotropic in spin space: therefore 
the spin-orbit coupling produces a spin diffusion current 
that tends to make the distribution isotropic. Our re- 
sults imply that the spin polarization s v created by the 
electronic current is significantly affected by spin-orbit 
coupling with the impurities, and differs from the stan- 
dard Edelstein^ and Aronov and Lyanda-Geller— pre- 
diction(see section 3.4). 

III. THE BOLTZMANN EQUATION 

In the following we will discuss the spin Hall effect in 
the context of the Boltzmann equation. We start with 
the standard Boltzmann equation in the absence of spin- 
orbit coupling, for which of course no spin Hall effect 
exists. In the subsequent steps intrinsic and extrinsic 
spin-orbit coupling will be included. In the extrinsic case 
we will concentrate on the side jump, that is we will ig- 
nore skew scattering terms. It will be demonstrated that 
in the presence of spin-orbit coupling a charge current 



generates a transverse spin current, i.e. a spin Hall effect 
appears. In the Rashba model, however, a voltage in- 
duced spin polarization also appears, whose contribution 
to the spin current exactly compensates the directly gen- 
erated spin Hall current. In other words, the very well 
known cancellation of the spin Hall effect in a (purely) 
Rashba 2DEG is in fact more general, and holds even in 
the presence of extrinsic spin-orbit coupling mechanisms. 

We relegate technical details to the appendices, which 
can be skipped in a first reading. 

A. No spin-orbit coupling 

To fix the notation we begin by considering the Boltz- 
mann equation when a — and Ao = 0, i.e. in the absence 
of spin-orbit interaction. We have 

dtf P + - • V x /p - eE • Vp/p 

m 

= -2w 2 ]>>( e P ~ £ p')(/p - fp'l (15) 
p' 

with e p = p 2 /2m and / p = /(p,x, t) is the distribution 
function. In the following / p will be assumed to have a 
spin structure in order to describe charge and spin dy- 
namics. We follow the notation of Ref. [2J and introduce 
the spin dependent density and current density as 

* — ' * — ' m 
p p 

The particle density is given by n — Tr(p), the spin den- 
sity by s a = Tr(cr a/ o)/2, a = x,y,z, and the particle 
and spin currents are defined analogously. Integrating 
the Boltzmann equation over the momentum p gives the 
continuity equation 

dtp + V x • J = 0, (17) 

i.e. the continuity equations for the charge and spin com- 
ponents, 

3 t n + V x -j=0, 9 4 s a + V x -j a = 0. (18) 

B. Intrinsic spin-orbit coupling 

Next, we consider the Rashba spin-orbit interaction in 
the absence of extrinsic mechanisms, i.e. we set Aq = 
but allow a ^ 0. This case has been discussed several 
times in the literature. For a review one may see Ref. l26l . 
Here we follow the treatment of Ref.[23|. Using the SU (2) 
formulation the Boltzmann equation becomes 

(a t + £.v x + i { F.v p ,})/p 

= ^w-o'E^p ~ £ p')(/p - /p')- (19) 
p' 
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Compared to the standard Boltzmann equation there 
are two modifications. First, in Eq. (fH?|) . the spacial 
derivative becomes the covariant derivative expected for 
a SU (2) vector potential A 

V x = V x + i[A,...], (20) 

where [,] denotes a commutator. Note that here the 
SU (2) vector potential includes only the Rashba contri- 
bution, i.e. only the first term on the r.h.s. ofEq.©. The 
second modification concerns the force F which appears 
in the anticommutator {, }, since in addition to the elec- 
tric field the force contains now a spin-dependent 5*1/(2) 
Lorentz-like force, 

F = — eE — — x B, B k = -e kln i [A u A n ] . (21) 
m l 

In the Rashba model the only non vanishing component 
of the "SU(2) magnetic field" is B z = — 2(ma) 2 a z . From 



The Boltzmann equation is modified and becomes 



the Boltzmann equation one finds immediately the con- 
tinuity equation 

dtp + V x ■ J = 0, (22) 

which has the same form as Eq. (TT7|) but for the spa- 
cial derivative, here replaced by the covariant one. The 
covariant derivative does not change the equation for 
the density, but there appear extra terms in the one for 
the spin density: the explicit result looks exactly like 
Eq. ((Til). 



C. Extrinsic spin-orbit coupling 

In the next step we consider the so-called side jump ef- 
fect due to extrinsic spin-orbit coupling, while we neglect 
the intrinsic contribution, Ag ^ 0, a = 0. 



dtf P 



ii/ p + 2 Wo 2 J2 - V)f 5 {(P ' - P) X /p ' } 
p' 

, r \2 1 r 



- eE • Vp/p 



/p-/p'-f i{(- e )E.(p'-p)x<r, 9/p ' 



9e r 



(23) 



The extra terms due to spin-orbit coupling appearing in 
Eq. (|23|) have been derived within the framework of a 
microscopic, Keldysh-based approach along the lines of 
Refs.Qjl and [I?]. More explicitly, the extra terms are ob- 
tained from the spin-orbit correction to the self-energy 
(Born approximation), compare Refs. [l6l and I27I and Ap- 
pendix [B] The term on the l.h.s. of the equation appear- 
ing under the spacial derivative corresponds to a cor- 
rection to the current operator, so that the continuity 
equation (|17p remains valid provided the expression for 
the current density is modified as follows 



J 



E 

p 



1 



£ /p + T ^{Px^/pl 

m 4 2r 



(24) 



We notice that both new terms contain the side jump 
shift Ax = — (Aq/4)Ap x <t, where Ap = p' p. Such an 
expression of the side jump can be derived by analyzing 
the effect of the spin-orbit interaction on the scattering 
trajectory of a wave packed. 

Spin-orbit coupling is assumed to be small, and in the 
following we will calculate the spin Hall current to low- 
est order in the parameter Ag , as usual in linear response 
to an homogeneous applied electric field. Spin-orbit cou- 



pling generates a spin Hall current in two ways: (1) due 
to the modified expression for the current operator and 
(2) due to the spin-orbit correction to the distribution 
function / p . 

The spin-orbit correction to the current operator - the 
second term on the r.h.s. of Eq. (p4|) - yields the spin 
current 



[3 



Ql(l) _ 



Ag m . 

8 T 



(25) 



where jb is the particle current in the absence of spin- 
orbit coupling. Expressing the latter in terms of the 
Drude conductivity, jf, — —(enT/m)E},, the result reads 



(1) _ 1 enA| 

Ui J — 2 4 



;bEb 



1 



<iab ^ Csj* E b- 



(26) 



Eq. (121)1) represents the first half of the side jump contri- 
bution to the spin Hall effect. 

The second half is found when one considers the first 
term on the r.h.s. of Eq. (|24j) and computes the spin- 
orbit correction to the distribution function / p . To this 
end, we multiply Eq. (|23l) by Pi/m, integrate over the 
momentum and consider the a-th spin component. 
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The last term on the l.h.s. of Eq. (|23[) would give a term quadratic in the electric field, therefore only the r.h.s. 
contributes and one is left with 

= -1-kthVq <H e P ~ £ p') 5— " 27rr?,iU ^ <5(e p - ep')-j-e bfcm (-eEb) K . (27) 

Tit £ ~r *-'^'T> > TTt Zi 

p,p' p,p' 



There follows 



N-°l(2) -TrV 7 " P f - lmX h p- e l rT sH W 



(28) 

By combining Eqs. (|26p and ([25} the total side jump 
contribution to the spin Hall conductivity reads a^f 1 = 



D. Intrinsic and extrinsic spin orbit coupling 

We are now ready to consider the interplay of intrin- 
sic and extrinsic spin-orbit coupling, i.e. both a and A§ 
are nonzero. Since we assume weak spin-orbit coupling 
one could be tempted to write a Boltzmann equation in 
which the spin-orbit terms of (|23[) and (|19[) are simply 
added. Such an equation was obtained in Refs. and 
I27I however it does not correctly capture the side jump 
contribution the the spin Hall conductivity and leads to 
the wrong result <r sH = (l/2)cr*j . Technically, what is 
missing in the derivation of the kinetic equation in [l6| and 
I27I (and also in the diagrammatic calculation of the spin 
Hall conductivity of Il3f) is the modification to the self- 
energy when intrinsic and extrinsic spin-orbit coupling 
are both present. An advantage of the present 517(2) ap- 
proach is that such a correction appears quite naturally, 
as discussed in Appendix [Bj The bottom line of includ- 
ing such a modification is that the correct Boltzmann 
equation is obtained by replacing in Eq. (|23p the spacial 
derivative with the covariant one, the force with the sum 
of the standard electric force and the SU(2) force, and by 
allowing for a correction to the collision integral due to 
the interplay of the two types of spin-orbit interaction. In 
this way one obtains again the continuity equation (|22[) 
- this time with J given by Eq. (|24p - which implies the 
vanishing of the spin Hall conductivity. 

The diffusive regime allows to analyze how this hap- 
pens in an analytically explicit way. Let us consider first, 
for simplicit y's sake, the pure Rashba case and set Ao = 0. 
We refer to [23| for details of the calculations. The cur- 
rent density can be written as the sum of a diffusion, 
drift, and Hall current J = Jdift + Jdrift + Jhaii with 



-L>V x p, Jdrift = crE, 



Jhall = 77- B X J. 

Ira 



(29) 

Here D = ijVpT is the diffusion constant, and a the elec- 
trical conductivity. Notice that the "Hall conductivity" 
is taken from the classical formula for the weak magnetic 



field regime of^ = <jdlj c t, where ao is the Drude con- 
ductivity and w c is the cyclotron frequency associated 
with the "SU(2) magnetic field". In the homogeneous 
case the explicit result for the spin current jy is 



3* = 2maDs v + a&E s 



(30) 



The first term has its origin in the diffusion current, while 
the second term appears due to the Hall current. One has 



r sH 



2t 



57T T DP 



(31) 



where top is the Dyakonov-Perel spin relaxation time, 
1/tdp = (2ma) 2 D. ' 

Since we now know that, due to the continuity equa- 
tion, the spin current j* must be zero, we can conclude 
that there must exist a spin polarization in y-direction, 
such that the diffusion current associated with it cancels 
the Hall current af^E x . Eq. ([30]) then implies 

1 ^ sH TP 



2maD 
-eaN rE x , 



(32) 



which is the well known result for the voltage induced 
spin polarization in the Rashba mode l 24 ' 25 . This may 
well be considered a new derivation of the Edelstein ef- 
fect. However, as we now show, the result for s y changes 
significantly when extrinsic contributions to the spin Hall 
current are taken into account. 

When extrinsic spin-orbit interaction is present too, 
the spin Hall current in the diffusive limit becomes 



3l = 2maD S y + 



)Ex, 



(33) 



so both the intrinsic and the extrinsic mechanisms drive 
a spin Hall current, but this is in the end compensated by 
the spin polarization s y . Notice that the vanishing spin 
Hall conductivity implies that the limit a — > does not 
reproduce the result obtained by setting a = from the 
outset. This paradox has a simple solution as pointed out 
infill: extrinsic spin-orbit coupling introduces the Elliott- 
Yafet spin relaxation mechanism, which is order Ag and 
thus does not appear in our equations so far. Elliott- 
Yafet spin relaxation is included by replacing Eq. (flj by 



d tS v + L s v = -2maj* 



(34) 



with 1 /t s the Elliott- Yafet spin relaxation rate given by 
(see, for instance, Ref. l27t ) 



Ts ~ T { 2 



(35) 
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From Eqs. (|33| and (|34|) one obtains 



1 



t s /tdp 
2ma 



sH 
int 



1/t s + 1/tdp 



(36) 



Eq. (|3l?)) shows that the Elliott- Yafet spin relaxation in- 
deed cures the non-analyticities in the limit a — > 0. 

Finally, including the skew scattering contribution to 
the spin Hall current can be done by simply adding it to 
Eq. (02]) to give 



31 



2maD S y + « + ai" + )E X . (37) 



For completeness in Appendix [C] we discuss the physi- 
cal origin of the skew scattering contribution and its mi- 
croscopic evaluation. When using Eq. (f37)) . rather than 
Eq. (|3"3")1 . into Eq. (jMJ), one obtains the same expressions 
as in Eq. (|3"6"]l with the replacement — > <r|j + . 
Observe that the expression (|36l) for s y is significantly 
different from the standard expressio n 24 ' 25 . The stan- 
dard expression is recovered only in the limit in which 
1/t s — > and the extrinsic contribution to a is ne- 
glected. 



to unify the description of intrinsic and extrinsic effects. 
Though for clarity's sake, and to establish a direct con- 
tact with previous works, we have focused on the Rashba 
model and thus on 5/7(2) fields, our approach can be 
used for any Hamiltonian with a linear-in-momentum in- 
trinsic spin-orbit term - possibly slowly varying in space 
and time, too. In such a unified picture the spin-charge 
coupled dynamics is naturally described by a set of sim- 
ple continuity equations. In particular, this allowed us 
to clarify the issue of the side jump contributions to the 
spin Hall effect in the presence of Rashba spin-orbit in- 
teraction. Furthermore, we have derived a new formula, 
Eq. ((36)) . for the transverse spin polarization associated 
with the electric current. This is significantly different 
from the Edelstein and Aronov and Lyanda-Geller for- 
mula and is amenable to experimental verification. 

We acknowledge useful discussions with Dimitrie Cul- 
cer. This work has been supported by the Deutsche 
Forschungsgemeinschaft through SPP1285 and by the 
French Agence Nationale de la Recherche through grant 
No. ANR-08-BLAN-0030-02. GV acknowledges support 
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IV. SUMMARY 



Appendix A: Equation of motion for the spin 
density — some technical details 



The interplay between intrinsic and extrinsic spin-orbit 
coupling effects in 2DEGs is an experimentally relevant 
issue whose theoretical description has proven rather del- 
icate and riddled with puzzles. Moreover, the latter has 
suffered from the lack of a physically transparent formu- 
lation. 

We have shown how formulating spin-orbit coupling 
in terms of non-Abelian gauge fields can be exploited 

I 



In the main text we showed that the equation of motion 
for the spin density can be written as 



d t s a + ddt = e abc A\jt + e abc 5A\5jl 



(Al) 



We mentioned that the correlations between the vector 
potential and the current vanish without giving any de- 
tails. Some will be given in this appendix. 



We begin by noting that with the introduction of the covariant momentum defined by 

n, = -\di + x -A\o h 

the paramagnetic and diamagnetic currents can be combined together to give for the total current 

«y?(x) = ^Tr]T m + U\),a c } < C t cp,)^'-")* 



(A2) 



(A3) 



where the density matrix (c^c p <) is a matrix in spin space. 
The fluctuating vector potential can be written as 



^(x)=i^e ii6 ^ gi 5U(x)e^ 



(A4) 



We then need to evaluate the average over impurities 

. An 



<Mf(x)tf?(x) Tr^{(H x n') 6 ,a c }JU(p-p')(4 V ). 
p p' 



(A5) 
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In obtaining the above equation, we have used th e fact that the standard white-noise disorder correlations, 
SV(x)SV(x') = mvlS(x — x'), require 5V(q)SV(q') = 7JiU 2 5(q + q'), so that qj in Eq. (|A4j) must be equal to 
Pj — p' j , which in turn is equal to IT, — U'j . We may then use the identity 

e ijb qj {n 4 + U' h <j c } = 2 {(n x W) b , a c } 

and obtain Eq. (|A5|) . 

To perform the average in Eq. (IA5|) we write the density matrix in terms of the Keldysh Green function to read 

<L4*(x)tf?(x) = -in iUo 2 ^Tr£ / ^ {(II x II') b , a c } (G£(e)G*(e) + G^G^e)) . (A6) 

p,p' 



Then, by observing that to first order in pp • A/(mei?) 

^nifflG^(e)=0, (A7) 
p 

the cancellation mentioned in the main text follows. 



Appendix B: Side jump corrections in the 
Boltzmann equation 

We determine the Boltzmann equation for a system 
with intrinsic and extrinsic spin-orbit coupling following 
the procedure of H^. The starting point is the Dyson 
equation for the (nonequilibrium) Green function from 
which the equation of motion 



9( + -'Vx-e-'E9 e 



>,v p .} 



G = -i[t,G] 
(Bl) 

for the "locally covariant" Green function G is obtained. 
t and x are the center-of-mass coordinates of the two- 
point Green function G and e, p are the Fourier trans- 
formed relative coordinates. Furthermore G and E have 
a structure in Keldysh space, i.e. they have retarded, ad- 
vanced and a Keldysh components. The locally covariant 
Green function, G, is related to the conventional Green 
function G by the transformation 

G(p) = G(p)-i{A-V p ,G(p)}, (B2) 

where A is the SU (2) vector potential. The distribution 
function appearing in the Boltzmann equation is given 
by the Keldysh component integrated over the energy, 



/p = 



1 



1 



(B3) 



so that by integrating (|B1|) over the energy one obtains 
the Boltzmann equation for the distribution function / p . 

Generally, the r.h.s. of (|B1I) is responsible for the scat- 
tering kernel in the Boltzmann equation. For example, 
within the Born approximation and for spin-independent 
scattering the self-energy is 



S = n 4 « 2 ^G(p) 



(B4) 



and leads to the r.h.s. of Eq. f(T9[) . 

The side jump corrections are found from the Born 
self-energy to first order in the (extrinsic) spin-orbit cou- 
plin g, A n- The expression for the self-energy is given in 
Ref. [27] for a = 0. In the presence of Rashba spin-orbit 
coupling the covariant expression, obtained by using the 
transformation (lB2l . reads E 1 = Ei + El + E\ with 



(B5) 



A 2 

p' 



An 



[pa<?bP' c ,G(p',X.) 

Ej(p,x) - ^n iUo 2 e a , c ^{a b ^,(V x ) a G(p',x)},(B6) 

E£(p, x) = -^n^e a6c £(V x i„ \ P ,n h . Gip'. x S f 137 ) 

p' 



where (V x ) a is the (spin) a-componcnt of the covari- 
ant derivative defined in Eq. ([20|) in the main text. No- 
tice that the covariant derivatives appearing in the self- 
energies Yi\ and E b act once inside and once outside the 
anticommutator. This is because the vector potential 
A a does not commute with the Pauli matrices o\>. When 
a = 0, the covariant derivative reduces to the simple spa- 
tial derivative and Eqs. (|B5IIB7[) reduce to the expressions 
reported in Ref. [2l]. In this case both and Ej can be 
written as divergence terms and taken to the l.h.s. of 
the kinetic equation (jBll) . giving rise to the spin-orbit 
corrections appearing on the l.h.s. of Eq. (l23l) . 

In the presence of an electric field the space derivative 
has to be replaced by 



V x -> V x - eE<9 £ , 



(B8) 



which gives rise to the electric field-dependent term on 
the r.h.s. of Eq.®. 

Finally, we note that although E^ is the leading term, 
it is not relevant for the side jump contribution the 
to spin Hall effect and has therefore been ignored in 
Eq. (|23[) . Indeed, this term is responsible for the so-called 
swapping of the spin currents 2 ^. 

When Rashba spin-orbit coupling is present, insert- 
ing the self-energies (|B5IIB7I) into the kinetic equation 
(|Bip leads to a Boltzmann equation where the side jump 
corrections to the current appear under the covariant 
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derivative, confirming the covariant form of the conti- 
nuity equation discussed in Section 2. In addition, on 
the r.h.s. of the Boltzmann equation the additional term 

p' 

appears. This can be interpreted as a correction to the 
collision integral arising from the interplay of extrinsic 
and intrinsic spin-orbit coupling. Such term, which is 
also manifestly covariant, is important for the consistency 
of the kinetic equation. However it can be easily seen 
that, due to the symmetry between p and p', it vanishes 
after summation over p and hence does not modify the 
continuity equation. 



Appendix C: Skew scattering 

For completeness we now discuss skew scattering. Let 
us assume a system with time reversal symmetry, but 
spin dependent scattering due to spin-orbit interaction. 
The scattering amplitude reads 

S = A + p x p' • <tB, (CI) 

where p and p' are unit vectors in the direction of the 
momentum before and after the scattering event. By con- 
sidering the scattering probability proportional to l^ 2 , 
one obtains three contributions given by \A\ 2 , \B\ 2 and 
2 IZe (AB*)p x p' • <r. Whereas the first two contribu- 
tions are spin independent and give the total scattering 
time, the third one represents the so-called skew scat- 
tering term according to which electrons with opposite 
spin are scattered in opposite directions. To lowest order 
in perturbation theory or Born approximation one has 
A — vq and B = — [(XqP^/AJvq, having adopted a point- 
like impurity potential V(x) = vqS^x). Clearly, since A 
and B are out of phase, there is no skew scattering effect 
to this order. For it to appear to first order in the spin- 
orbit coupling constant A 2 ,, A has to be evaluated beyond 
the Born approximation. The scattering problem can be 
cast in terms of the Lippman-Schwinger equation 

V>(x) = e ikx + J dx' G(x - x')V r (x')V'(x'), (C2) 

where G(x) is the retarded Green function at fixed en- 
ergy. From (IC2[) we get 

= v G(x),A {1) = v ; 
?/> (2) = v 2 G(0)G(x),A^ = ugG(0). (C3) 



Notice that only the imaginary part of A^ is needed. 
By recalling that Im G(0) = —ttN , the spin-orbit inde- 
pendent scattering amplitude reads 

A = v (l- mN a v ) . (C4) 

The skew scattering contribution will then follow by in- 
serting the modified scattering amplitude (|C4j) into the 
collision integral of the Boltzmann equation. The same 
result can, of course, be obtained in quantum field theory 
using the Green function technique. The latter becomes 
necessary when one wants to consider skew scattering 
in the presence of Rashba spin-orbit interaction. To this 
end one has to consider the electron self-energy at least to 
third order in the scattering potential vq. The expression 
for this was obtained in Ref. [271 for the case a = and 
gives rise to three terms depending on the position of the 
insertion of — (Aq/4)ct x VV^(x) • p into a third order di- 
agram. When Rashba spin-orbit coupling is present one 
has to consider the covariant self-energy, as done for the 
side jump contribution in appendix [Bj It may be shown 
that to leading order, this is done simply by replacing the 
Green function G with its covariant expression G. Hence 
the self-energy responsible for the skew scattering reads 
E ss = E" s + Ef + Ef , where 

E* s = -mvl^f J2 6 (p°) ^fa) Pfr x p ■ <r (C5) 

Pa,Pt 

Ef s = -imv^ J2 3(P«) P« x Pb ■ a- G( Pb ) (C6) 

Pa,P6 

Ef = -imvl^ J2 PXPaV G( Pa ) G( Pb ). (C7) 

Pa,Pt, 

Here the standard impurity average has been performed, 
so that translational invariance is recovered and the self- 
energy depends on the external momentum p only. Since 
we are considering the effect to first order in Aq, the co- 
variant Green functions entering Eqs. ([C5IIC7p are spin in- 
dependent and isotropic in momentum space. As a result 
the retarded and advanced components of the above self- 
energies vanish, while the Keldysh component survives 
only for E^ s and E^f. Their joint contribution, after re- 
calling that ^2 G R (p) = —ittNo, leads to an extra term 
on the r.h.s. of the Boltzmann equation 

A 2 

- 2im l v 2 (voirNo)-2-^26(e p - e p >) {p' x p • a, f p ,} . 

(C8) 

Its contribution to the spin current is obtained as usual 
by multiplying Eq. (IC8|) by Pi/m, integrating over mo- 
mentum and projecting on the a-th spin component. 
Again, (|C8j) is already order Aq, so the distribution func- 
tion f p > is the one in the absence of spin-orbit coupling. 
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One has 

-2irniv£{v irN )— > 6(e p - e p /)— e jlk p ; pi 2/", = e iaj - {v TrN ) E j} C9) 

p'.p 

and the expression appearing in Eq. ([3]) follows at once. 
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